Quantum criticality as a resource for quantum estimation 
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We address quantum critical systems as a resource in quantum estimation and derive the ultimate quantum 
limits to the precision of any estimator of the coupling parameters. In particular, if L denotes the size of a 
system and A is the relevant coupling parameters driving a quantum phase transition, we show that a precision 
improvement of order 1/L may be achieved in the estimation of A at the critical point compared to the non- 
critical case. We show that analogue results hold for temperature estimation in classical phase transitions. 
Results are illustrated by means of a specific example involving a fermion tight-binding model with pair creation 
(BCS model). 
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I. INTRODUCTION 



It is often the case that a quantity of interest is not directly 
accessible, either in principle or due to experimental impedi- 
ments. In these situations one should resort to indirect mea- 
surements, inferring the value of the quantity of interest by 
inspecting a given probe. This is basically a parameter es- 
timation problem whose solution may be found using tools 
from classical estimation theory |Q]] or, when quantum sys- 
tems are involved, from its quantum counterpart [2]. Indeed, 
quantum estimation theory has been successfully applied to 
find optimal measurements and, in turn, to evaluate the cor- 
responding lower bounds on precision, for the estimation of 
parameters imposed by both unitary and nonunitary transfor- 
mations. These include single-mode phase-shift flflELdis- 
placement @], squee zing 10, IH as well as depolarizing 131 or 
amplitude-damping channels in finite-dimensional sys- 
tem, lossy channel in infinite-dimensional ones |[Tll[T2Tl . and 
the position of a single photon 111 311 . Here we focus on the 
estimation of the parameters driving the dynamics of an inter- 
acting many -body quantum system i.e., the coupling constants 
defining the system's Hamiltonian and temperature. It is a 
generic fact that when those are changed the system is driven 
into different phases. When this happens at zero temperature 
one says that a quantum phase transition (QPT) as occurred 
0]. Different phases of a system are, by definition, char- 
acterized by radically different physical properties, e. g. the 
expectation value of some distinguished observable (order pa- 
rameter). This in turn implies that the corresponding quantum 
states have to be statistically distinguishable more effectively 
than states belonging to the same phase. It is a non trivial 
result that the degree of statistical distinguishability is quanti- 
fied by the Hilbert-space (pure states) distance 11511 or by the 
density operator distance (mixed states) ifl^l . This amount to 
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say that the Hilbert-space geometry is an information-space 
geometry ifnl [H, O . One is then naturally led to consider 
the distance functions between infinitesimally close quantum 
states obtained by an infinitesimal change of the parameters 
defining the system's Hamiltonian. Boundaries between dif- 
ferent phases are then tentatively identified with the set of 
points where this small change of parameters gives rise to a 
major change in the distance i.e., major enhancement of the 
statistical distinguishability. This is precisely the strategy ad- 
vocated in 12(X. l2ll I22L I23L l24l 12511 in the so-called metric (or 
fidelity) approach to critical phenomena. 

The purpose of this paper is to establish the following re- 
sult: quantum criticality represent a resource for the quantum 
estimation of Hamiltonian parameters. Indeed, by exploiting 
the geometrical theory of quantum estimation, we will show 
that the accuracy of the estimation of coupling constants and 
field strengths at the critical points is greatly enhanced with 
respect to the non-critical ones. We will also discuss ultimate 
limits imposed by quantum mechanics to this scheme of pa- 
rameter estimation. These results are, under several points of 
view, analogue of those showing that entanglement is a useful 
metrological resource 1261,1271 12811 . In particular, if L denotes 
the size of a system and A is the relevant coupling parame- 
ters driving a quantum phase transition, we will show that a 
relative improvement of order L may be achieved in the esti- 
mation of A at the critical point. 

The paper is structured as follows: In Section[II]we provide 
a brief introduction to the geometrical theory of quantum esti- 
mation whereas in Section [Til] we review the metric approach 
to quantum criticality. In Section [IV] we show how quantum 
critical systems represent a resource for quantum estimation 
and derive general results about ultimate quantum limits to 
precision. In Section [V] we illustrate properties of the opti- 
mal measurements, also for systems at finite temperature, and 
the connection with the optimal measurement for state dis- 
crimination. In Section[VT]we illustrate our general results by 
means of a specific example involving a fermion tight-binding 
model with pair creation. Section I VIII closes the paper with 
some concluding remarks. 
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II. GEOMETRY OF QUANTUM ESTIMATION THEORY 

The solution of a parameter estimation problem amounts to 
find an estimator, i.e a mapping A = X(x\,X2, ■■■) from the 
set x °f measurement outcomes into the space of parameters. 
Optimal estimators in classical estimation theory are those sat- 
urating the Cramer-Rao inequality, 

V X [X] >F-\\) 

which poses a lower bound on the mean square error 
V\ [X]jk — E\ [(A — A) j (A — X)k] in terms of the Fisher infor- 
mation 

F(A) = / dX(x) P CX\X)[d x logpCX\X)} 2 . 



square error of any quantum estimator is bounded by the in- 
verse of the quantum Fisher information. In formula 



v x [x] > n-\x) 



(i) 



QCR is an ultimate bound: it does depend on the geometrical 
structure of the quantum statistical model and does not depend 
on the measurement. Notice that due to noncommutativity of 
quantum mechanics the bound may be not attainable, as it is 
the case for several multiparameter models. When one has at 
disposal M identical copies of the state qx the QCR bound 
reads Vx[X] > 1/MH(X), which is easily derived upon ex- 
ploiting the additivit y of the quantum Fisher information. A 
relevant remark lTl6l l29ll is that the SLD itself represents an 
optimal measurement and the corresponding Fisher informa- 
tion is equal to the QFI. 



Of course for unbiased estimators, as those we will deal 
with, the mean square error is equal to the covariance matrix 
V x [X] jk = Ex[XjX k } - Ex[Xj]Ex[Xk] . 

When quantum systems are involved any estimation prob- 
lem may be stated by considering a family of quantum states 
q(X) which are defined on a given Hilbert space Ti and la- 
beled by a parameter A living on a d-dimensional manifold 
M., with the mapping A i— » g(X) providing a coordinate sys- 
tem. This is sometimes referred to as a quantum statistical 
model. In turn, a quantum estimator A for A is a selfadjoint 
operator, which describe a quantum measurement followed by 
any classical data processing performed on the outcomes. As 
in the classical case, the goal of a quantum inference process 
is to find the optimal estimator. The ultimate precision at- 
tainable by quantum measurements in inferring the value of 
a parameter or a set of parameters is expressed by the quan- 
tum Cramer-Rao (QCR) theorem lfl6l l29ll which sets a lower 
bound for the mean square error of the quantum Fisher infor- 
mation. The QCR bound is independent on the measurement 
and very much based on the geometrical structure of the set 
of the involved quantum states. The symmetric logarithmic 
derivative C(X) (SLD) is implicitly defined as the (set) Her- 
mitian operator(s) satisfying the equation(s) 

MA) = i{e(A)r M (A) + £ / ,(A) e (A)} , 

where := d/d^, (fi = 1, . . . , d). From the e above equa- 
tion, when Qj + Qk > 0, one obtains 

(<Pj\£nW\<Pk) = 1{ ( Pj\ d >*Q( > >)\<Pk) I {Qj + Qk), 

where we have used the spectral resolution g(X) — 
^2k Qk\<fk){<Pk\- The quantum Fisher information H(X) 
(QFI) is a matrix defined as follows 



T~tfj,v{X) = Tr 



e {X)l{jC ti {X)C u {X)+£ v (X)C ti {X)} 



The QFI is symmetric, real and positive semidefinite, i.e. rep- 
resents a metric for the manifold underlying the quantum sta- 
tistical model U3,[Hl- The QCR theorem states that the mean 



HI. GEOMETRY OF QUANTUM PHASE TRANSITIONS 

Now we turn to illustrate the achievements of the metric 
approach to QPTs that are relevant to the problem investi- 
gated in this paper. The crucial point is that this approach, 
being based on the state-space geometry is universally appli- 
cable to any statistical system and does not require any pre- 
liminary understanding of the structure of the different phases 
e.g., symmetry breaking patterns, order parameters. In prin- 
ciple not even the system's Hamiltonian has to be known as 
long as the relevant states are (see e.g., the analysis of matrix- 
product state QPTs [24(]). Of course the main conceptual as 
well as technical obstacle one has to overcome in order to get 
this simple strategy at work is provided by the fact that the 
relevant phenomena are intrinsically thermodynamical limit 
ones. More technically, one considers the set of Gibbs ther- 
mal states 00(A) := Z- l e-P H ( x \ (Z := Tr[e- 0H ^}) asso- 
ciated with a parametric family of Hamiltonians {H(X)}x^m- 
Physically the A's are to be thought of as coupling constants 
strengths and external fields defining the many-body Hamilto- 
nian H(X). The systems one is interested in are characterized 
by the fact that, in the thermodynamical limit, they feature 
a zero-temperature i.e., quantum, phase transitions (QPT) for 
critical values A c fl4ll . We now consider the Bures metric 
ds 2 B = y"V q av dX^dXu over the manifold of density matri- 
ces where [301 



It. 



(<Pj\d l j.Q\ ( Pk)('Pk\dvQ\'Pj) 
Qk + Qj 



(2) 



The key result of the metric (or fidelity) approach to QPTs is 
that the set of critical parameters can be identified and ana- 
lyzed in terms of the scaling and finite-size scaling behavior 
of the metric ©. More precisely the metric has the following 
properties: i)In the thermodynamical limit and in the neigh- 
borhood of the critical values A c the zero-temperature metric 
has the following scaling behavior ds 2 B ~ L d \X — A c |~ 9 , 
where L is the system size, d the spatial dimensionality, v 
the correlation length exponent (£ ~ |A — X c \~") and A g = 
2£ + d — 2 Ay. Here £ is the dynamical exponent and Ay is 
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the scaling dimension of the operator coupled to A. ii) At the 
critical points, or more generally in the critical region defined 
by L <C £, the finite-size scaling is as follows ds 2 B ~ L d + As . 
The main point is that for a wide class of QPTs A g can be 
greater than zero, whereas at the regular points the scaling is 
always extensive i.e., ds 2 B ~ L d . To be precise superextensive 
behavior requires that the perturbation be sufficiently relevant. 
The superextensive behavior gives rise, for L — > oo, to a peak 
(drop) of the metric (fidelity) that allows one to identify the 
boundaries between the different phases. Moreover it can be 
proven [31] that, for local Hamiltonians, superextensive be- 
havior of any of the metric elements is a sufficient condition 
for gaplessness i.e., criticality. On the other hand criticality its 
not a sufficient condition for such a super-extensive scaling. 
There are indeed QPTs driven by local operators not suffi- 
ciently relevant (renormalization group sense) where no peak 
in the metric is observed at the critical point [31]. iii) When 
the temperature is turned on one can still see the signatures of 
quantum criticality. This is done by studying the scaling, as a 
function of the temperature, of the elements of the metric (O 
or fidelity [32]. In particular when the temperature is small 
but bigger than the system's energy gap one has ds 2 B ~ T~P, 
((3 > 0). When one sits, in the parameter space, at the criti- 
cal point this result can be extended all the way down to zero 
temperature giving rise to a divergent behavior that matches 
L\X - \c\-" As - Remarkably, crossovers between semiclas- 
sical and quantum critical regions in the (T, A) plane can be 
identified by studying the largest eigenvalue of the metric or 
its curvature lO . 



ample if A = /xq — A', then A' can be effectively estimated 
around A' = po — A c (po is assumed to be exactly known). 

In order to quantitatively asses the improvement in the es- 
timation accuracy, let us focus on the single parameter case 
A £ R. In this case the QCR reads V\ [A] > (4gn) _1 . From 
the results summarized in the previous section one easily sees 
that i) In the neighborhood of critical point and at zero tem- 
perature the optimal covariance V\ [A] scales like | A — A c \ I/As . 
This is a remarkable fact for at least two reasons. On the one 
hand, it means that the covariance itself scales as the param- 
eter and thus divergence of the QFI are allowed IU2I1 . On the 
other hand, for those QPTs such that A 9 > 0, the covariance 
can be then pushed all the way down to zero by getting closer 
and closer to the quantum critical point, ii) The covariance of 
A may achieve the limit L~ a where a = d in all regular i.e., 
non-critical points, while at the critical ones one can achieve 
a > d. For example in the class of one-dimensional systems 
studied in 11221 l23l 12411 one has that the estimation accuracy 
goes from order L^ 1 in the regular points to L~ 2 at the critical 
points, iii) At the critical point in the parameter space and for 
finite temperature T the covariance of the optimal estimator 
scales like T°. The exponent (3 is related to the A 9 above and 
for a class QPTs is greater than zero 113 311 . In this case by ap- 
proaching zero temperature accuracy grows unboundedly, i.e 
QFI diverges II 1211 . In contrast, at the regular points accuracy 
remains finite even for T — > (and finite system's size). In 
the quasi-free fermionic case mentioned above one has (3=1. 



A. Finite-size corrections 



IV. CRITICALITY AS A RESOURCE 

To the aims of this paper it is crucial to notice that the QFI 
is proportional to the Bures metric (0. Indeed, by evaluat- 
ing the trace defining the QFI in the eigenbasis of g(X) one 
readily finds fl^l <^„ = jH^. This simple remark, along 
with the results of the metric approach to criticality summa- 
rized in the previous section, immediately lead to the main 
conclusion of this paper: the estimation of a physical quantity 
driving a quantum phase transitions is dramatically enhanced 
at the quantum critical point. It is important to notice that 
the Hamiltonian dependence on the quantity to be estimated 
can be even an indirect one. More precisely, suppose that H 
depends on a set of coupling constants A and those in turn 
depend on the unknown (to be estimated) quantities A' i.e., 
A = /(A'), through a known function /; we assume also that 
/ is smooth and that its derivatives are bounded and system's 
size independent. From the tensor nature of the Bures metric 
i.e., g'^ v = gapidXa/dX'^idXp/dXl), one has that the QFI 
associated to the A' has the same dependence on the system's 
size of the QFI associated to the A's and the same divergen- 
cies in the thermodynamical limit. From the QCR bound (fl~|i 
then it follows that if one is able to engineer a system such 
that the coupling constants A defining its Hamiltonian (featur- 
ing QPTs) depends on the unknown quantities A' in the way 
outlined above then the A"s can be estimated with a greater 
efficiency in the points corresponding, to the QPTs. For ex- 



In view of practical applications involving realistic sam- 
ples, we now consider corrections arising from the finite sys- 
tem size. As we have seen, when L — » oo, the maximum of 
the QFI is located at the critical point A c . Instead, for finite L, 
the location of the maximum is shifted by an amount which 
goes to zero as L — > oo. To obtain an estimate of the shift one 
must include off-scaling contribution. The previous formu- 
lae for the scaling of ds B , TL in the off- and quasi-critical re- 
gions can be combined in a single equation valid in a broader 
regime: 



g = L A ^(z) + L^D(X) + 



-Ao-< 



C (A) + smaller terms. 



(3) 



Here z is the scaling variable z = (L/^) 1 ^, e > 0, and <j> 
is a scaling function satisfying 4>{z) ~ z~ l/Aa when z — > oo 
(the off-critical region), whereas we must have </>(0) ^ to 
comply with the behavior in the quasi-critical region (z — > 0). 
The existence of such a scaling function is a consequence of 
the scaling hypothesis. Instead the functions C(A), D(X) are 
analytic around A c and are responsible for corrections to scal- 
ing. The maximum of H(X) defines a pseudo-critical point A^ 
whose location, for large L is given by 

X l A c « ~^jjL-^ ^L- 2 ^ ZjL- 2 '^. (4) 
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One can say that the shift exponent of the QFI is, \jv, 2/v, 
or larger depending on the form of the functions above. A 
similar shift of the location of the quantum critical point may 
be observed when the temperature is turned on. 



V. OPTIMAL MEASUREMENTS 

Our results establish the quantum limits to the degree of 
accuracy in estimating the coupling constant of a locally- 
interacting many-body system. It is crucial to stress again that 
the QCR bound is, in this case, attainable with the correspond- 
ing SLD representing an observable to be measured in order 
to achieve the optimal estimation. Of course the SLD may be, 
in general a very complex observable which itself depends on 
the ground or thermal state of the system. For a family of pure 
states g(X) = \ip(X))(ip(X)\, the SLD is easily derived as 

£„(A) = 20„ e (A) = 2(1^(^1 + \d^)W) ■ 

Using first order perturbation theory, one finds (r^d^ip) — 
(E n — Eq)^ 1 (n\d^H 0} where now the |n)'s {En's) are the 
eigenvectors (eigenvalues) of H (A) being |0) the ground state. 
Therefore, we may write the SLD as 

C fl (X)=2[(P d tl HG(E ) + h.c] 

where P = \0)(0\, G{E) = Pi[H(X) - E]~ 1 P 1 , being P x = 
I — Pq. In the one parameter case the non- vanishing eigen- 
values of the SLD are given by ±2y/(drp\d^) — \(^p\dip)\ 2 — 
±2dsB /d\ from which the unknown parameter A can be esti- 
mated. The corresponding QFI is given by B25I1 



?V(A)=4^ 



n>0 



(O|a M ff|n)(n|0„ff|O) 
(E n - E ) 2 



(5) 



From this expression one sees that the origin of the divergence 
of QFI at the critical point is the vanishing (in the thermody- 
namical limit) of one of the E n — Eq factors. 



A. Finite temperature 

Now we would like to show that also classical i.e., tem- 
perature driven, phase transitions provide in principle a valu- 
able resource for estimation theory. Consider the metric 
induced on the (inverse) temperature axis by the mapping 
f3 — * g((3). In M25I1 and 113 311 it has been shown that ds 2 ~ 
df3 2 ({H 2 )(3 - (H) 2 ) = df3 2 p- 2 c v ((3), where c v (^denotes 
the specific heat at the inverse temperature (3. In B25I1 it was 
noticed that this relation suggests a neat and deep interplay be- 
tween Hilbert space geometry and thermodynamics. Here we 
stress that also quantum estimation gets involved. Indeed, fol- 
lowing the same lines used in the QPTs case, one easily real- 
izes that when the specific heat shows an anomalous increase 
then the same happens to the QFI associated to the parameter 
f3. From this fact stems that at the classical phase transitions 
with diverging specific heat one can estimate temperature with 
arbitrarily high accuracy. Conversely, if the specific heat is 



bounded from above the estimation accuracy of the tempera- 
ture is bounded from below. Again, in analogy with the QPT 
case discussed above, if can be made dependent, in some 
known fashion, on some other parameter A', then this latter 
can be estimated with better accuracy at the phase transition. 
In other words, we have a quantitative statement of the intu- 
itive expectation about the fact that thermometers are more 
precise at the points where changes of states of matter occur. 



B. Quantum discrimination 

In our analysis the problem of interest is that of estimat- 
ing the value of a parameter or a set of parameters. On the 
other hand, when one deals with discrimination of parameters 
rather than estimation, the relevant metric dsn n R can be de- 
rived from the so-called quantum Chernoff bound [34], which 
arises in the problem of discriminating two quantum states in 
the setting of asymptotically many copies. In view of the in- 
equality \ds 2 B < dsg CB < ds 2 B all the results obtained here 
can be generalized to dsg CB . This means that quantum criti- 
cality is a resource for quantum discrimination as well. 



C. Applications: two-stage adaptive measurements 

A question may arise on how our results may be exploited 
in practice, being the form of the SLD, which maximizes the 
Fisher information, typically dependent on the true, unknown, 
state of the quantum system. The apparent loophole in the 
argument is closed by noticing that one can still achieve the 
same rate of distinguishability by a two-stage adaptive mea- 
surement procedure [35j]. Roughly speaking the estimation 
scheme goes as follows: one starts by performing a generic, 
perhaps suboptimal, measurement on a vanishing fraction of 
the copies of the system and obtains a preliminary estima- 
tion. Then one measures the remaining copies, taking the 
the preliminary estimation as the true value. This guaran- 
tees that the Fisher information obtained by the second se- 
ries of measurements approaches the QFI as the number of 
measurements goes to infinity, and that the resulting estima- 
tor saturate the QCR bound. Notice that the achievability of 
the QCR bound is ensured when a single parameter has to be 
estimated 13611 . though the actual implementation of the opti- 
mal measurement, which is general not unique ll37ll . may be 
more or less challenging depending on the specific features 
of the system under investigation. We also notice that any 
maximum-likelihood (ML) IBal estimator, defined as the es- 
timator A = X(x\,X2, ■■•) maximizing the likelihood function 
C(X) = Ylk p( x k\X), M being the number of measurements 
and p(x\X) the conditional probability density of the measured 
quantity, is consistent, i.e. it converges in probability to the 
true value and asymptotically (M 3> 1) efficient, i.e. it sat- 
urates the Cramer-Rao bound in the limit of many measure- 
ments, thus achieving the ultimate bound in precision. 
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VI. QUANTUM ESTIMATION IN THE BCS MODEL 

In order to appreciate the critical enhancement of precision 
in a specific physical situation we consider a fermionic tight- 
binding model with pair creation i.e. the BCS-like model de- 
fined by the Hamiltonian 



where 



b z (d) = l e- ikd b z 



kesz 



and analogously for b y (d). For example, for L = 4 the general 
formula reduces to 



H, 



J 



i=l 



4<h + i + j44 , 



h.c.) - 2hJ2c\ Cl 



(6) 



where L denotes the system size and periodic boundary con- 
ditions have been used. Upon considering the thermal or the 
ground state of the above Hamiltonian we have a uniparamet- 
ric statistical model where J is the parameter to be estimated, 
h is the external (tunable) field and the anisotropy 7 is a fixed 
quantity. For 7^0 this model undergoes a quantum phase 
transition of Ising type at h = J (because of symmetries 
there is an analogous critical point at h = —J). Precisely 
around the Ising transition point point we will investigate the 
enhancement in precision offered by criticality. 

Moving to Fourier space the Hamiltonian may be rewritten 

as 



keBZ 



1)+ {^kc\cl k +h. 



keBZ 



= ^[-e k Tj> + A k T»] 



(7) 



k>0 



where e k = — J cos (k) — h, A k = — J7 sin (fc), the Brillouin 
zone BZ ranges from — ir to tt and the momenta are of the 
form k = 2-nn/L, with integer n. The expression © for 
Hj follows from the observation that in the subspace spanned 



by ||0), <4 c -fc|0)} the operators n k + n_ fc 
h.c. 



1 



and 



( 



lc k c -k 



- j» 



k , represent a set of Pauli operators tL 

(and they are zero in the complementary space). 

In the above quasi-spin formulation the Hamiltonian in 
Eq. (0 has a block form, and so is the thermal state 
q,j = Z^ 1 exp{— {3Hj}. Being each block essentially two- 
dimensional, the SLD of the model can be obtained by com- 
puting the SLD in each block, with the aid of formula (18) of 
039(1 ■ One then arrives at 



£(J) = J2KTk+btr 



y 

k ' 



fc>0 



where, at T = 



h Afcefc 



h -JAl 



(8) 



where A k = yf? k 



Going back to real space we have 
1 

2 



£(J) = ^L6*(0)-5>t&* (l-j) Cj 



1 J 



1,] 



04+h.c. 



(9) 



C(J) 



" (^ 2 +7 2 J 2 ) 3 / 2 
J 7 



J 7 



c|c3 + c|c4 + h.c. 



J-fN 
2 

+ h 



(10) 



t J 



C[C2 ~ C\C 



1°4 



which represents a collective measurement on the system. 
TV = J2j c ] c j is the total number operator. 

In general the coefficients b z (d) and b y (d) decay exponen- 
tially at normal points of the phase diagram and thus the SLD 
is a local operator. On the other hand, the decay is only alge- 
braic in the critical region \h — J\L < 1 and this corresponds 
to a SLD given by a collective measurement. Indeed, at the 
Ising transition one can show that b v (d) ~ d~ l . For b z (d) the 
integral of Fourier coefficient does not converge for large L, 
so that one has to keep the sum with L finite. The sum is well 
approximated by b z (d) ~ 2.9 |£ — 5I Notice that for 

7 = the SLD at zero temperature is identically zero since 
for any finite L changing J only results in a level crossing. 

When the temperature is turned on, one can draw similar 
conclusions. The local character of the optimal measurement 
is enhanced at regular points VL whereas in the critical regime 
nonlocal measurements are needed. More specifically, the co- 
efficients b z ' y (d) decay exponentially at regular points and al- 
gebraically in the region \h — J\ < T. 

The ultimate precision is determined by the QFI Hj, which, 
for such quasi-free Fermi system reduces to (see 12511 ) 



k>0 



k>0 



h 2 "f 2 sin (ky 



(11) 



where = arctan (e k /A k ). Upon introducing the scal- 
ing variable z = L (h — J) one is interested in the behavior 
around the Ising transition, i.e. for z <C 1 where one can 
observe superextensive behavior. Indeed, upon expanding for 
small z and using Euler-Maclaurin formula, one finally gets 



H(J) 



L 2 



L 



24J 2 7 2 

„2 



Z 

L 

fl 

L 2 



r 



2tt 2 J 2 7 2 
1) L 2 



12J 3 7 4 
L 4 



720JV 



+ 0(L) 
0(L 2 ) 



+ O (z 3 ) 



(12) 



From this expression one can explicitly see the improvement 
of precision due to the underlying quantum critical behavior of 
the system. Besides, one can locate the pseudocritical point, 
defined as value of the field leading to the maximum ofH(J). 
Differentiating the above formula with respect to h one ob- 
tains 

h* L = J + 30J ( 7 2 - 1) L-' 2 + O (L- 3 ) , 
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and use this information to achieve the quantum Cramer-Rao 
bound also for finite-size systems. 

With the notations of section lTV Al we note we can write the 
scaling form of Tt( J)/L as in Eq. (f3]l with A g = v = e = 1 
and 



4>{z) 



l 



24JV 720JV 



0(z 3 ) 



h 2 - 1) r 2 

d (h) = ^2J^r ( h -J) + [(h - J) 

C(h)^^^0(h-J) . 



(13) 
(14) 
(15) 



Having 4>'(0) = the shift exponent turns out to be 2/v = 2 
At finite temperature and in the thermodynamic al limit the 
QFI has been calculated previously for quasi-free Fermi sys- 
tem lf33ll . In the present case we obtain 



H(J) _P 2 



dk 



(J + h cos (k)Y 



8tt J cosri ((3A k /2) 



(16) 



1 

2^ 



77 cosh(/3A fc ) - 1 fry Bin (A) 



cosh(/3A fc ) 



■dk. (17) 



In this case we verified numerically that the maximum of 
H(J) always occurs at h = J where it has a cusp. Then, 
since when the temperature goes to zero the first integral van- 
ishes we evaluate the second term at the critical point. For 
h = J the dispersion Afc is linear around k = n, which gives 
the dominant contribution to the integral. One then obtains 



2C 

H(J) = - 2 



TT 2 T\Jj\ 



+ O (T°) 



C = 0.915 being the Catalan constant. From the above for- 
mula one can again appreciate the enhancement of the bound 



to precision occurring when T goes to zero in the critical 
regime. 



VII. CONCLUSIONS AND OUTLOOKS 

In conclusion, upon bringing together results from the ge- 
ometric theory of quantum estimation and the geometric the- 
ory of quantum phase transition we have quantitatively shown 
that phase transitions represent a resource for the estimation of 
Hamiltonian parameters as well as of temperature. To this aim 
we used the quantum Cramer-Rao bound and the equivalence 
of the notion of quantum Fisher metric and that of quantum 
(ground or thermal) state metric. We have also found an ex- 
plicit form of the observable achieving the ultimate precision. 
A specific example involving a fermionic tight-binding model 
with pair creation has been presented in order to illustrate the 
critical enhancement and the properties of the optimal mea- 
surement. The improvement in estimation tasks brought about 
by quantum criticality is reminiscent of the one associated to 
quantum entanglement in computational as well as metrologi- 
cal tasks. The analysis reported in this paper makes an impor- 
tant point of principle and establishes the ultimate quantum 
limits to the precision with which one can estimate coupling 
constants characterizing a quantum Hamiltonian. Since the 
Hamiltonian completely characterize the quantum dynamics 
one can say that our results shed light on the ultimate lim- 
its imposed by quantum theory to the observer capability of 
knowing a quantum dynamics. Remarkably the boundaries 
between different phases of quantum matter are where these 
limits gets looser and a deeper knowledge can be achieved. 

MGAP would like to acknowledge useful discussions with 
Alex Monras and Paulina Marian. 
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